Scalable simultaneous multi-qubit readout with 99.99% single-shot fidelity 
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We describe single-shot readout of a trapped-ion multi-qubit register using space and time-resolved 
camera detection. For a single qubit we measure 0.9(3) x 10~ 4 readout error in 400 /xs exposure 
time, limited by the qubit 's decay lifetime. For a four-qubit register (a "qunybble") we measure an 
additional error of only 0.1(1) x 10 -4 per qubit, despite the presence of 4% optical cross-talk between 
neighbouring qubits. A study of the cross-talk indicates that the method would scale with negligible 
loss of fidelity to ~ 10000 qubits at a density < 1 qubit / /im 2 , with a readout time ~ 1 /is/qubit. 

PACS numbers: 03.67.-a, 37.10.Ty, 42.50.Dv, 42.30.Va 



Single-shot qubit readout is an essential requirement 
for quantum computing, in particular for implement- 
ing quantum error-correction [l| or measurement-based 
quantum computing 0. Fault-tolerant quantum com- 
puting requires readout fidelities >99.9% if a large in- 
crease in the number of physical qubits per logical qubit 
is to be avoided [3j. Two major issues in scaling up 
from single-qubit to many-qubit systems are the need 
to perform operations on many qubits in parallel [l| , and 
to overcome errors due to "cross-talk" between different 
qubits. In this paper we describe simultaneous readout 
on a multi-qubit register and show that the quantum 
state can be accurately inferred even in the presence of 
significant cross-talk. Since the latter is short-range, the 
method is scalable to much larger numbers of qubits. 

For any physical implementation, readout error is usu- 
ally limited by the time available to integrate a detec- 
tor signal (to improve the signal-to-noise ratio) before 
the qubit decays. In recent work on a two-electron spin 
qubit in a double quantum dot, for example, the trade- 
off was between electronic measurement noise and the 
34 /us qubit lifetime Cross-talk in a multi-qubit sys- 
tem could increase the measurement noise and/or de- 
crease the qubit lifetime. In trapped-ion or neutral-atom 
approaches to quantum computing, qubit readout is usu- 
ally achieved by driving a fluorescing cycling transition 
which involves one of the qubit states but not the other, 
and measuring whether or not the atom fluoresces ||. 
Here the noise is predominantly photon-counting shot 
noise and the qubit lifetime is limited by spontaneous 
decay or off-resonant excitation, which can transfer the 
qubit from the "dark" state to the "bright" state or vice 
versa. Previous work reported measurement fidelities of 
90%-98% for two- and three-qubit systems 

The highest fidelity quantum logic gate demonstrated 
to date [9J] used a trapped-ion optical qubit stored in 
the (4S 1/2 , 3D 5/2 ) levels of 40 Ca+, fig. QJi. For the 
same qubit, we have already described single-shot read- 
out with 99.991(1)% fidelity for a single ion, using time- 
resolved detection of photons counted with a photomul- 
tiplier tube (PMT) [lOj. Here, we extend this to simul- 
taneous readout of multiple qubits using an electron- 
multiplying charge-coupled device (EM-CCD) camera, 
which allows both space- and time-resolved detection. 



We first measure the readout error for a single qubit 
and find that, as with the PMT, at any given photon 
collection efficiency this is limited by the decay lifetime 
t = 1168(7) ms of the metastable D 5/2 level [ll|- Then we 
measure the readout error for a four-ion "qunybble" [l2| , 
where we eliminate the effect of the D 5/2 decay by post- 
selecting data, to find the residual error due to cross-talk 
between the imperfectly imaged ions. 

The linear ion trap apparatus (Tlj and the methods 
for preparing and detecting the bright Si /2 and dark D 5/2 
qubit states are similar to those in [Tpj . with the PMT 
replaced by the camera. Fluorescence at 397 nm from the 
ion(s) is collected by a compound objective of numerical 
aperture sin a = 0.25 and imaged onto the camera with 
net efficiency 1.0(1)%. The objective is not optimized for 
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FIG. 1: (a) 40 Ca+ levels, (b) Image cross-talk affecting ion 2 
in a string of four ions (inset: 50 x 10 pixel acquisition area, 
at scale 2.6 /im/pixel). Filled circles show the integrated fluo- 
rescence from ion 2; the decreasing-brightness pixel order and 
data were obtained by averaging images in which only ion 2 
was fluorescing. For the other three ions (hollow symbols) the 
same pixels were taken in the same order, but the data were 
taken from images in which only that ion was fluorescing. Up 
to ~ 100 pixels, the region-of-interest (ROI) is roughly cir- 
cular, centred on ion 2, with a diameter in the object plane 
given by the upper abscissa. The lines show calculated per- 
formance for a diffraction-limited imaging system of the same 
numerical aperture, with ion spacings ten times smaller than 
in the image, on a 10 x magnified upper abscissa. 
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397 nm. resulting in performance well below the diffrac- 
tion limit (fig-lb)- The nominal quantum efficiency (QE) 
of the camera at 397 nm is 48%, approximately twice that 
of the PMT. However, the stochastic gain process in the 
camera's EM register leads to an excess noise factor of V2 
above shot noise when operated at high gain, equivalent 
to an effective reduction in QE by a factor 2 for photon 
counting applications [13j. Background counts are dom- 
inated by the camera's clock-induced charge (CIC) read- 
out noise rather than by scattered light, but the finite 
qubit lifetime is still the principal limitation to readout 
fidelity. Cosmic rays cause fewer errors for the CCD than 
the PMT [l(J, due to its spatial discrimination. 

We define the average readout error per qubit as e = 
■^( e B + 6 d)i where cb is the fraction of experiments in 
which an ion prepared in the bright state was detected 
to be dark, and similarly for eo- In the single-qubit ex- 
periments described below, there is a small population 
preparation error for the dark state [10j . resulting in a 
contribution toeofO.I3(l)xfCF 4 ; values of e given below 
are after subtraction of this quantity. In the four-qubit 
experiment the state preparation error is expected to be 
negligible because post-selection is used. 

The readout fidelity for a single ion was measured by 
repeatedly preparing and measuring each qubit state, 
and comparing the measurement result with the known 
state preparation, for a total of 103744 trials. The 
method is similar to that described in [r|, except that 
the time-resolved PMT detection bin is replaced by a 
single 400 /is camera exposure. Camera images were 
analysed both by simple pixel-count thresholding and 
by maximum likelihood techniques. We label the pixels 
i = 1 ... N in order of decreasing average brightness (de- 
termined from a long-exposure image, fig. [2ji) and select 
a region-of-interest (ROI) by varying N. 

For the threshold analysis method, we histogram the 
total pixel counts n i where n% is the photon count 

for pixel i, separately for the bright and dark state prepa- 
rations, and choose the count threshold which optimizes 
discrimination between the histograms (fig. [2Jd, inset). 
The variation of e with ROI size N is shown in fig. [2b; it 
decreases with TV as the histogram overlap falls, reaches 
a minimum e = 0.9(3) x 10~ 4 (with e B = 0) for N = 28, 
then eventually rises when the ROI is too large due to in- 
creasing background counts from pixels far from the ion 
which contain negligible information. 

The simple threshold method does not take into ac- 
count the known spatial distribution of the ion fluores- 
cence over the ROI. In the "spatial" maximum likelihood 
analysis, we consider the likelihood pb = Il^Ii Bii n i) 
that the set of pixel counts {n-i} came from a bright ion, 

and compare this to the likelihood po = Yli=i^i{ n i) 
that {n-i} arose from a dark ion. Here Bi(rii) is the prob- 
ability of observing rii counts on pixel i for a bright ion; 
similarly for Di. Since the probability distribution Bi 
is different for each pixel, the information contained in 
the spatial distribution of the fluorescence is used. The 
distributions Bi and Di can be obtained from indepen- 



SPACE-RESOLVED TIME — RESOLVED 




1 10 100 500 2 4 6 8 10 12 14 16 18 

no. pixels AT (in brightness order) no. of 200jiis exposures 



FIG. 2: Single-ion readout results for space- and time-resolved 
detection. Analysis methods: (T)hreshold, (M)aximum like- 
lihood, (A)daptive maximum likelihood, (a) Integrated flu- 
orescence signal and long-exposure image, (b) Single 400 fis 
exposure: readout error e versus ROI size N (inset: bright and 
dark photon count histograms for N = 30, with threshold 6). 
(c) Time-resolved detection using multiple 200 fis exposures 
with, above, example exposures 1-5 for bright and dark state 
preparation (inset: minimum e versus N for method M). 



dent control data, or by fitting the distributions from the 
experiment itself; either approach gave similar results. 
Maximum likelihood analysis results are also shown in 
fig-Hb; the minimum value of e = 1-0(4) x 10 -4 is achieved 
with an ROI of only N = 10 pixels, demonstrating that 
the maximum likelihood technique is more efficient than 
the threshold method. However, it does not lead to a 
lower e as this is limited by the D 5/2 decay. 

An adaptive detection method (analogous to that de- 
scribed in (7l.[Tfj|) can reduce further the average number 
of pixels N required to determine the qubit state. The 
estimated probability of measurement error can be in- 
ferred from pb and pu- For reliable readout, we only 
need use sufficient pixels for the absolute log-likelihood 
ratio R = \ \ti{pb/pd)\ to be above some chosen confi- 
dence level. Taking the pixels in brightness order en- 
sures most efficient use of the information. This method 
reaches the same e in N = 2.9 pixels, fig. EJj. 

We also investigated time-resolved detection [lj] for a 
single qubit. In this experiment, the single 400 /xs cam- 
era exposure was replaced by a series of 18 exposures, 
each of duration t s = 200 /xs, to investigate whether time- 
resolved detection could give an improvement in readout 
fidelity by detecting any of the D 5/2 decays which lim- 
ited the single-exposure fidelity. The method is similar 
to that used for PMT detection in [r|, except that the 
CIC noise associated with reading each exposure from 
the CCD means that it is disadvantageous to make each 
exposure too short. The CCD itself was used to store all 
18 exposures, allowing a rapid sequence of images with 
minimal dead time (6 fis) between each. Alternate bright 
and dark state preparations gave a total of 67892 trials 
(« 1.2 million exposures), which were analysed by simple 
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thresholding and by maximum likelihood analysis taking 
into account the effect of D 5/2 decay, with results shown 
in fig. [2k. For the "spatio-temporal" maximum likelihood 



(<0 



(b) approx. circular ROI diameter at ions (/^m) 



to m 10 mo a>'- 



analysis we compared the likelihoods pb — Y\jLi PBj and 
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which allows for the probability (t s /r) of spontaneous 
decay in exposure j' of M (see [Hi). The probability 
PBj that exposure j came from a bright ion is given by 
Pbj = TliLi -Bji(riji) where Bji(riji) gives the probabil- 
ity of measuring riji counts on pixel i of exposure j for 
a bright ion; similarly for PDj- The minimum error is 
e = 1.1(4) x 10~ 4 , for M > 4 exposures, no lower than 
the single-exposure experiment. A (temporally) adaptive 
analysis is also shown in fig. [2}:, and reaches the same e in 
an average of M = 1.1 exposures (230 ^s); however the 
overhead in CCD read time for multiple exposures means 
the adaptive technique offers little practical speed-up. 

We now describe an experiment to measure the readout 
error ex for the four-ion qunybble, independently of the 
D 5/2 decay error. The cross-talk due to imperfect imaging 
in our system is illustrated in fig. [TJb, where the data 
were taken from the experiment to be described below. 
For an ROI centred on ion 2 with diameter equal to the 
mean ion spacing of 14 /im, the nearest neighbour ions 
each contribute signals w 4.0% of the signal from ion 2 
itself, while the next-nearest neighbour contributes 0.9%. 
Cross-talk at this level could clearly impact the readout 
error at the 10 -4 level achieved for a single ion. 

We are not able to prepare all 16 possible qunybble 
states {0000, 1111} deterministically with an error 
~ 10~ 5 as in the single-qubit experiments above. Instead, 
in each trial we prepare a random qunybble state, and 
sandwich the readout test measurement between "pre" 
and "post" measurements. We then select for analysis 
only those trials in which the "pre" and "post" measure- 
ments agree: these measurements give the state that was 
prepared, to which the "test" measurement is compared. 
A dark— ^bright decay from the D 5/2 level during the trial 
will cause the "pre" and "post" measurements to disagree 
for the ion which decayed. By increasing the exposure 
time of the "pre" and "post" measurements only, we in- 
crease the accuracy of this state post-selection procedure 
at the expense of rejecting more trials. 

A single trial consists of 6 identical camera exposures 
each of 400 fis duration (fig. IS 1 ), with the CCD being read 
out between each. An initial "check" exposure is taken 
with all ions fluorescing. Then a short (12 fxs) pulse of 
393 nm light shelves the ions in D 5/2 with a probabil- 
ity 0.46 per ion, to prepare a random qunybble state 
(with each state approximately equally probable, inset 
hg- 13b) ■ The next two exposures are added to give the 
"pre" measurement, the fourth exposure constitutes the 
actual readout "test" measurement, and the last two ex- 
posures are added for the "post" measurement. Finally 
an 854 nm laser pulse returns all ions to the bright state. 
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FIG. 3: Four-ion readout experiment, (a) Experimental se- 
quence, showing a typical trial; maximum likelihood anal- 
ysis applied to the 400 /j,s "test" exposure gives the qun- 
ybble state as {1110} with estimated probability of error 
~ Sfc=i e ~ Rk = 1-4 x 10 -11 . (b) Readout error per qubit av- 
eraged over trials of all 16 qunybble states {0000, . . . , 1111}, 
whose distribution is shown in the inset. The error due to 
D 5 / 2 decay is excluded as described in the text. Analysis 
methods: (T)hreshold; (M)aximum likelihood ignoring neigh- 
bours; (MN) iterative maximum likelihood method taking 
into account nearest neighbours, with (solid line) a Monte- 
Carlo simulation of 10 9 trials. Adjacent ROIs overlap when 
their diameters exceed the 14 )im mean ion spacing. 



This sequence was repeated 30100 times, giving 120400 
qubit trials. We first identify those trials where the "pre" 
and "post" measurements disagree, and discard them, us- 
ing thresholding of pixel counts summed over rectangular 
ROIs. Two thresholds were used to improve our confi- 
dence in the retained trials: we demand either that the 
"pre" and "post" measurements are both above the up- 
per threshold, or that both are below the lower threshold. 
The separation between the upper and lower thresholds 
was such that 95% of trials were retained [ill ]. 

Having thus selected a set of trials with known qubit 
states, we proceed to analyse the "test" exposures. We 
first apply a simple threshold analysis, treating each ion 
independently. The pixel counts n k i for each ion k are 
summed over independent brightness-order ROIs, and an 
independent threshold is chosen for each ion. The min- 
imum readout error (per qubit) is ex = 6.8(8) x 10~ 4 
(fig. [3)3) and ex rises rapidly above TV sa 25 pixels, where 
the ROI diameters exceed the mean ion spacing. 

The spatial maximum likelihood analysis should be 
more immune to the effects of image cross-talk since it re- 
quires fewer pixels and gives less weight to pixels further 
from the ion of interest. For each ion k we compare the 
likelihoods p B = Ilili B ki (n kl ) and p D = \[f =l D ki (n k i) 
that the set of counts {n k i} came from a bright or dark 
ion, where B k i{n k i) is the probability of observing n k i 
counts on ion fc's ith brightest pixel when that ion is 
bright, and similarly for D k i(n k i). The probability dis- 
tributions B k i and D k i are obtained from the "pre" and 
"post" exposures only and are thus independent of the 
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"test" exposures. Results are shown in fig. |3Jd and reach 
an error ex = 0.7(3) x 1CP 4 per qubit, an order of magni- 
tude below simple thresholding. As expected, the error is 
less sensitive to the ROI overlapping neighbouring ions. 

Nevertheless, we can improve the analysis by taking 
into account the effect of neighbouring ions explicitly. 
There is a significant effect on the dark distribution Dk 
depending on whether the neighbour k' is bright or dark; 
indeed for a pixel far from ion k but close to ion k' 
the shape of the distribution will be dominated by the 
state of k! even though the pixel may still contain use- 
ful information on the state of A:. To extract this infor- 
mation, we now allow independent probability distribu- 
tions B v ki(jiki) and D v ki{nki) for each neighbour state 
{v} = {00, Of , f 0, f 1 }. These distributions are again 
found using the "pre" and "post" exposures only, where 
we know the state of both ion k and its neighbours {f }. 
However, when we use these distributions to analyse a 
particular "test" exposure we do not assume any knowl- 
edge of the neighbour states: instead we proceed iter- 
atively, making an initial guess for the qunybble state 
of {0000}, measuring the state of each qubit k with 
{is} = {00}, recalculating using the neighbour states 
from the first iteration to give a second measurement, 
and so on until the inferred qunybble state is stable. 
Results from this iterative maximum likelihood analysis 
are also plotted in fig. \3jp: for ROI size N > 60 pix- 
els the readout error per qubit reaches an average level 
of e x = O.f (f) x 10~ 4 and is even zero (in 1 14606 tri- 
als) for certain N. There is now no penalty for allowing 
the ROI to overlap neighbouring ions. The agreement 
with the simulation in fig.[3}D indicates that systematic ef- 
fects (e.g., correlated errors common to "pre", "post" and 
"test" exposures) are below the statistical uncertainty. 

Finally we can search for effects due to non-nearest 
neighbours by generalizing the preceding analysis to 
neighbour states {v} = {000, . . . , 111}, which allows for 
all permutations of the three neighbour ions to any ion 
k. We find no reduction in ex or N and conclude that 
any effect is at or below the ~ 10" 5 level. 

The single-ion experiment achieved a readout fidelity 
of 99.991(3)% with a single 400 [is camera exposure, lim- 
ited by the photon collection efficiency and the lifetime 
of the D B/2 qubit state. The four-ion qunybble experi- 
ment demonstrated that, when the effects of the D 5/2 de- 
cay are excluded, discrimination accuracy of 99.999(1)% 



can be achieved with the same exposure time, even in 
the presence of significant image cross-talk from neigh- 
bouring ions. Together these two experiments imply that 
99.99% net fidelity can be attained with the camera for 
simultaneous multiple qubit readout. Since it was found 
sufficient to include only the effects of nearest neighbors, 
this result should extend to one-dimensional qubit regis- 
ters of arbitrary length, given comparable imaging qual- 
ity over the entire ion string. As the error due to neigh- 
bour effects is small (ex <C e), the fidelity should be 
similar even for a 2D array of ions (or neutral atoms in 
a long- wavelength optical lattice [l5[ , given the same flu- 
orescence and dark-B-bright rates). This 512 x 512 pixel 
camera could be used to measure a 2D array of at least 
fOOOO qubits simultaneously; the exposure time would 
then be negligible compared with the CCD readout time 
(~ f /is per qubit at the fastest available readout speed of 
0.1 /xs/pixel). A CCD which allowed parallel pixel read- 
out would offer a significant speed-up. Image processing 
time is < O.f /is/pixel, so the main classical computing 
resource cost is storage (~ f kB/pixel) for the Bi and Di 
distributions. 

The cross-talk in our system was limited by imper- 
fections in the imaging system. To estimate what is 
theoretically possible, we plot on fig. QJ> the calculated 
image cross-talk for a diffraction-limited imaging sys- 
tem at A = 397 nm with the same numerical aperture 
sin a = 0.25, but with ion spacings reduced by a fac- 
tor of 10. The cross-talk from neighbouring ions was 
calculated by integrating the Airy pattern over off-axis 
circular ROIs, and is qualitatively similar to our mea- 
sured cross-talk. In this calculation the spacing between 
the ions, f.4 /jm, is less than the diameter of their Airy 
discs f.22A/ tana = 1.9 ^m. The implication is that 
readout fidelity at the 99.99% level is attainable for ions 
or atoms spaced close to the diffraction limit, with an 
imaging system possessing the same collection efficiency. 
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